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The inclusive production of charmonium mesons in B meson decay has been studied in a 20.3 fb~^ 
data set collected by the BABAR experiment operating at the T(45) resonance. Branching fractions 
have been measured for the inclusive production of the charmonium mesons J/^ , V'(2S'), Xci, and 
Xc2- The branching fractions are also presented as a function of the center-of-mass momentum of 
the mesons and of the helicity of the J/ip . 

PACS numbers: 13.25.Hw, 14.40.Gx 



I. INTRODUCTION 

Studies of the inclusive production of charmonium 
mesons in B decays provide insight into the physics of 
the underlying production mechanisms. Non-relativistic 
QCD (NRQCD) Q, which may provide an explanation 
^] for the unexpectedly large production of J/ip and 
'0(25') mesons observed in pp collisions is an exam- 
ple of such a mechanism. NRQCD calculations use phe- 
nomenological matrix elements that should be applicable 
to a variety of production processes 0], including such 
kinematically different regimes as hadron collisions at the 
Tevatron collider and e+e~ collisions at the T(4S'). 

This paper presents an analysis of J/tp , Xci, Xc2, and 
ipC^S) mesons produced in B decays at the T(45') res- 
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onance. J/ijj and il){2S) mesons are reconstructed in the 
e+e" and decay modes, and the Xci and Xc2 in 

the J/V'7 final states. We also reconstruct ip{2S) 
J/'i/;7r"*"7r~. 

The results include new measurements of previously 
observed decays 0, Q as well as momentum and helicity 
distributions not previously measured. 

II. THE BABAR DETECTOR AND THE PEP-II 
COLLIDER 

The BABAR detector is located at the PEP-II e+e" 
storage rings operating at the Stanford Linear Acceler- 
ator Center. At PEP-II, 9.0 GeV electrons coUide with 
3.1 GeV positrons to produce a center-of-mass energy of 
10.58 GeV, the mass of the T(4S') resonance. 

The BABAR detector is described elsewhere (7j; here 
we give only a brief overview. Surrounding the interac- 
tion point is a 5-layer double-sided silicon vertex tracker 
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(SVT), which provides precision spatial information for 
all charged particles, and also measures their energy loss 
(dE/dx). The SVT is the primary detection device for 
low momentum charged particles. Outside the SVT, a 
40-layer drift chamber (DCH) provides measurements of 
the transverse momenta pr of charged particles with re- 
spect to the beam direction. The resolution of the pr 
measurement for tracks with momenta above 1 GeV/c is 
parameterized as 

= 0.13 -^7% + 0.45%, (1) 

Pt 

where pr is measured in GeV/c. The drift chamber also 
measures dE/dx with a resolution of 7.5%. Beyond the 
outer radius of the DCH is a detector of internally re- 
flected Cherenkov radiation (DIRC), which is used pri- 
marily for charged hadron identification. The detector 
consists of quartz bars in which Cherenkov light is pro- 
duced as relativistic charged particles traverse the mate- 
rial. The light is internally reflected along the length of 
the bar into a water-filled stand-off box mounted on the 
rear of the detector. The Cherenkov rings expand in the 
stand-off box and are measured with an array of photo- 
multipher tubes mounted on its outer surface. A CsI(Tl) 
crystal electromagnetic calorimeter (EMC) is used to de- 
tect photons and neutral hadrons, as well as to identify 
electrons. The energy resolution of the calorimeter is pa- 
rameterized as: 

a{E) 2.3% , „^ 

^ = — ^ ® 1.9%, (2) 

where the energy E is measured in GeV. The EMC is 
surrounded by a superconducting solenoid that produces 
a 1.5-T magnetic field. The instrumented flux return 
(IFR) consists of multiple layers of resistive plate cham- 
bers (RPC) interleaved with the flux return iron. The 
IFR is used in the identification of muons and neutral 
hadrons. 

Data acquisition is triggered with a two-level system. 
The first level uses fast algorithms implemented in hard- 
ware that examine tracks in the DCH and energetic clus- 
ters in the EMC. The second level retains events in which 
the track candidates point back to the beam interaction 
region, or in which the EMC cluster candidates are cor- 
related in time with the rest of the event and exceed the 
energy of a minimum ionizing particle. Over 99.9% of 
BB events pass the second level trigger. A fraction of 
all events that pass the first level are passed through the 
second to allow monitoring of its performance. 



III. COORDINATE SYSTEM AND REFERENCE 
FRAMES 

We use a right-handed coordinate system with the z 
axis along the electron beam direction and the y axis up- 
wards, with origin at the nominal beam interaction point. 



The polar angle 9 is measured from the z axis and the az- 
imuthal angle from the x axis. Unless otherwise stated, 
kinematic quantities are calculated in the rest frame of 
the detector. The other reference frame we commonly 
use is the center of mass of the colliding electrons and 
positrons, which we call the center-of-mass frame. 



IV. EVENT SELECTION 

The data used in these analyses were collected between 
October 1999 and October 2000 and correspond to an 
integrated luminosity of 20.3 fb^^ taken on the T(45) 
and 2.6 fb~^ taken off-resonance at an energy 0.04 GeV 
lower than the peak, which is below the threshold for 
BB production and therefore includes only continuum 
processes. 

We use an equivalent luminosity of simulated data , 
including both BB and continuum, to study the effi- 
ciency of the analysis. 

We require events to satisfy criteria that are intended 
to have high efficiency for T(45') events while rejecting 
a significant fraction of continuum events and strongly 
suppressing beam gas events. The event must satisfy 
either the DCH components of both trigger levels or the 
EMC components of both, and have three or more high- 
quality tracks in the angular region with full tracking 
acceptance, 0.41 < 9 < 2.54. Reconstructed charged 
particles are considered high-quality tracks if they have 
at least 12 hits in the DCH, pT > 100 MeV/c, and a point 
of closest approach to the beam spot of < 3 cm in z and 
< 1.5 cm in xy. 

The ratio of the second to the zeroth Fox- Wolfram mo- 
ment i?2 measures how uniformly the energy in the 
event is distributed, distinguishing the more spherical 
BB events with small i?2 from the more jet-like contin- 
uum events at large i?2 (Fig.0. We require i?2 < 0.5. 

The primary event vertex is obtained from all charged 
particles with 0.41 < 9 < 2.54. Particles contributing 
a large a-re removed from the vertex and the process 
is iterated until stable. To reject events due to a beam 
particle striking the beam pipe or a residual gas molecule, 
we require the primary event vertex to be within 0.5 cm 
of the beam spot in xy and 6 cm in z. 

Finally, the event must include total visible energy 
(charged particles plus unassociated EMC clusters above 
30MeV) greater than 4.5 GeV. 

The efficiency for simulated T{4S) BB events is 
95.4%. More relevant is the ratio Ce of this efficiency 
to that for events containing the charmonium decay of 
interest. Ce is calculated for each of the final states us- 
ing simulated data. Inaccuracies in the simulation of the 
tracking systems produce an uncertainty of 1.1%, com- 
mon to all modes. 
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FIG. 1: i?2 distribution after all other selection criteria have 
been applied. Points are on-resonance data, histogram is off- 
resonance data scaled to the same luminosity. The vertical 
line denotes the requirement imposed on 7?2. 



A. Determination of Number of T(4S) Mesons 

The number of T{AS) events satisfying the above se- 
lection criteria (Nx) is obtained from the total number 
of events satisfying the criteria by subtracting the com- 
ponent due to the continuum. 



Nj 



Non-K-Roff - Mo 



(3) 



where 



Non is the number of events satisfying the criteria in the 
on-resonance data set; 

is the number of muon pairs in the on-resonance 
data set; 

Roff = ^off/Moff is the ratio of the number of events 
satisfying the hadronic selection criteria to the 
number of muon pairs in the off resonance data; 
and 

K = 1.0000 ± 0.0025 allows for differences in the ratios 
of continuum and muon-pair cross sections and ef- 
ficiencies between on-resonance and off-resonance 
data. 

The two highest-momentum tracks in muon pair events 
must both deposit less than 1 GeV in the EMC and satisfy 
I cos 6*1 < 0.7485 in the center-of-mass frame (to be within 
the region of full tracking efficiency) . The tracks must be 
within 10° of being back-to-back, and must have a com- 
bined mass greater than 7.5 GeV/c^. Since the number of 
muon pairs in the on and off resonance data sets appears 
only as a ratio in Equation O Nr depends only weakly 
on the details of the selection criteria. Varying these over 
reasonable ranges changes Nr by 0.5%, which we take as 
a systematic error. The on-resonance data contains 7.8 
times as many muon pairs as the off-resonance data. 
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FIG. 2: Mass distribution of J/tp candidates reconstructed in 
the (a) e"'"e~ and (b) fj,^ final states. The vertical dashed 
line in (b) marks the lower edge of the mass range used in the 
fi^ j-L~ final state. 



Changes to the trigger configuration caused Roff to 
vary from 4.89 to 4.94 during the period in which the data 
was collected. For the purposes of this calculation, the 
data is grouped into periods of compatible Roff- Com- 
bining the data in a single group changes Nr by 0.3%, 
which is taken as a systematic error. In principle, Roff 
could also vary due to beam gas backgrounds. However, 
the z distribution of the primary event vertex indicates 
that less than 0.1% of events are due to beam gas. 

We find Nr = (21.26 ± 0.17) x 10*^. The 0.8% un- 
certainty is systematic and includes the 0.5%, 0.3% and 
0.1% contributions from muon pairs, Roff, and beam 
gas described above. The largest component is from the 
0.25% uncertainty on k, which corresponds to a 0.6% 
uncertainty on Nr- 



V. J/i> PRODUCTION 

A. J/ip Reconstruction 

We reconstruct J/ip candidates in selected events us- 
ing the e'^e~ and ^i'^ final states. The leptons are 
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required to satisfy the track quality criteria listed earlier. 
Electron candidates are further required to have a dis- 
tance of closest approach to the beam line of less than 
0.25 cm to reject electrons produced by photon conver- 
sions. 

Both leptons are required to fall in the angular range 
0.410 < 9 < 2.409 rad (the overlap of the SVT and 
EMC coverage). Simulation indicates that (75.3 ±0.9)% 
of J/ip decays give both leptons in this region. As de- 
scribed in Sec. IIX Al the J/ip momentum distribution in 
the simulation 10] is slightly different from the observed 
distribution. This difference, and a small variation of ef- 
ficiency with momentum, produce the uncertainty in the 
efficiency. 

We obtain the efficiency for the leptons to satisfy the 
quality criteria by comparing the performance of the in- 
dependent SVT and DCH tracking systems in hadronic 
events. The corresponding uncertainty in the efficiency 
is 2.4% per J/ip . 

One particle in a J/tp e+e~ candidate must satisfy 
the "very tight" electron identification criteria described 
below. The other must satisfy the "tight" criteria dis- 
played in square brackets. 

• Difference between measured and expected energy 
loss in the DCH between —2a and +Aa, where a is 
the measurement error [—3(7 and +7 a]; 

• Ratio of energy measured in EMC to measured mo- 
mentum E/p in range 0.89 to 1.2 [0.75 to 1.3]; 

• Associated EMC cluster must include at least four 
crystals [same]; 

• Lateral energy distribution LAT of EMC clus- 
ter in range 0.1-0.6 [0.0-0.6]; 

• The A42 Zernike moment 12] of the EMC cluster 
< 0.42 [no requirement]; and 

• DIRC Cherenkov angle within 3a of expected value 
[no requirement]. 

LAT is a measure of the radial energy profile of the 
cluster, and is used to suppress clusters from electronic 
noise (very low LAT) or hadronic interactions (high 
LAT). A42 measures the azimuthal asymmetry of the 
cluster about its peak, distinguishing electromagnetic 
from hadronic showers. 

We reduce the impact of bremsstrahlung by combin- 
ing photons radiated by electron candidates with the 
track measured in the tracking system ( "bremsstrahlung- 
recovery" ) . Such photons must have EMC energy greater 
than 30 MeV and a polar angle within 35 mrad of the 
electron direction. The azimuthal angle of the photon 
must be within 50 mrad of the electron direction at the 
beamspot or be between this direction and azimuthal lo- 
cation of the electron shower in the EMC. 

For J/^p i-J.'^ IJ.^ candidates, one muon candidate 
must satisfy "tight" criteria, while the other satisfies 
"loose" criteria (shown in square brackets): 



• Energy in calorimeter between 0.05 and 0.40 GeV 
[< 0.50 GeV ]; 

• Number of IFR layers Nj pn > 2 [same] ; 

• Particle penetrates at least 2.2 interaction lengths 
A of detector material [2A]; 

• Measured penetration within ±0.8A of the value ex- 
pected for that momentum [±1.0 A]; 

• An average of less than 8 hit strips per IFR layer 
[same] ; 

• The RMS of the hits per IFR layer less than 4 
[same] ; 

• For candidates in the forward endcap, the number 
of hit IFR layers divided by the total number of 
layers between the first and last hit layers must be 
> 0.34 to reject beam background in the outermost 
layer [> 0.30]; 

• of the track fit in the IFR < 3 x Njpn [same]; 
and 

• of match between track from SVT and DCH and 

that found in the IFR < 5 x Njpji [same]. 

The particle identification efficiencies are obtained by 
comparing the yield of J/^ mesons applying various crite- 
ria to one or both tracks. The efficiency for J/ip e+e" 
satisfying the angular acceptance and track quality cri- 
teria is 90.5% with a systematic error of 1.8%. For 
J/ip it is 71.7% with a systematic error of 

1.4%. The systematic errors are somewhat conserva- 
tive, in that they include a component due to the J/4> 
statistics. Misidentification of hadrons as muons is higher 
than misidentification as electrons, producing J/^ back- 
ground levels that are approximately a factor of two 
higher. 

Finally, the J/ip candidate must have momentum in 
the center of mass p* < 2.0GeV/c. This requirement is 
fully efficient for a J/ip meson from B decay but rejects 
approximately 74% of those produced in the continuum 

More than one J/ip candidate may be found in an 
event. A second candidate is observed in 0.8% of events 
that include at least one. 



B. Extraction of Number of J/tp Mesons 

The mass of the J/tp candidate is obtained after con- 
straining the two leptons to a common vertex. In less 
than 1% of events, the vertex fit does not converge, and 
we instead use four-vector addition to obtain the can- 
didate mass. The mass resolution is poorer by approx- 
imately 1% for these events. Figure |21 shows the mass 
distribution of the selected candidates in the two lepton 
modes. 
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The number of J/ip mesons in the mass window used in 
the fit (2.6-3.3 GeV/c^ for the e+e" mode, 2.8-3.3 GeV/c^ 
for is determined by a binned hkehhood fit to the 

distribution. The background is represented by a third- 
order Chcbychev polynomial. The probability distribu- 
tion function (pdf ) for the /i^/i^ signal is the distribution 
from simulation convolved with a Gaussian distribution, 
with mean (allowing for a systematic shift between sim- 
ulation and data) and width (allowing for poorer resolu- 
tion) free to float in the fit. The additional smearing 
is required because the simulation underestimates the 
actual amount of material in the detector. The fit re- 
turns an offset of 3 MeV/c^ and an additional resolution of 
7.8MeV/c^. The total mass resolution in data is approx- 
imately 12MeV/c^. The simulation includes final state 
radiation and predicts that (2.7 ± 0.1)% of recon- 
structed J/ij} M^M^ candidates fall outside the mass 
range used in the fit. 

To include the impact of bremsstrahlung, the 
J/tp — !■ e+e~ signal pdf includes four components, corre- 
sponding to mesons where neither electron has undergone 
bremsstrahlung or at least one, and mesons for which the 
bremsstrahlung-recovery process has located a photon or 
not. The shapes are derived from simulated data, but 
the relative weights of three of the four components are 
allowed to float in the fit. The fraction of events that 
did not undergo bremsstrahlung but had a photon as- 
signed by the bremsstrahlung-recovery process is nomi- 
nally fixed to the value predicted by simulation, although 
it is varied to obtain a systematic error. 

An estimated (7.3 ± 0.8)% of reconstructed 
J/'if) e+e" candidates fall outside the mass range, 
which can be compared to the value of 6.1% if the 
relative pdf weights are fixed to the values predicted by 
simulation. 

We perform similar fits to the off-resonance data with 
all signal fit parameters fixed except for the number of 
mesons. The result is scaled by the ratio of on- to off- 
peak luminosity and subtracted to obtain the number of 
mesons attributable to B decay, which appears in Tabled 
as the net meson yield. 

The fitting procedure is validated and systematic er- 
rors on its results are obtained by comparing the gen- 
erated number of events with the fit number of events 
for many simulated mass distributions convolved with a 
Gaussian distribution. For the purposes of this test, we 
increase the statistics of the simulated data by relaxing 
the particle identification and track-quality requirements. 
We also test second and fourth-order Chebychev polyno- 
mials for the background pdf. We perform these tests for 
the Xci, Xc2, and 'ip{2S) mass distributions as well. 

We vary fit parameters that are fixed during the fit to 
obtain an additional systematic contribution. In the case 
of the J/ip , we vary the bremsstrahlung-recovery error 
rate from one-half to twice its nominal value. The sys- 
tematic errors on the fit yields are 0.7% for J/tp — > ijl^ ji'' 
and 4.1% for J/il) e+e~. 



Determination of the B 
Fraction 



J/tpX Branching 



We calculate values for the B — > J/ip X branching frac- 
tion using the e^e~ and n'^fJ-~ final states separately, 
then combine the two. The equation for the branching 
fraction is the same for both cases (and for the other 
mesons studied): 



2-Nr-ec-Bc 



Ce, 



(4) 



where 



is the net number of mesons in the mass fit range 
after continuum subtraction; 

ec is the efficiency for a meson to satisfy the selection 
criteria, including the requirement that the mass 
fall in the mass fit range; 

Be is world average 15] for the relevant secondary char- 
monium branching fraction. For the J/ip , this is for 
J/ip — > e+e~ or J/ip — * M'*'/^"; 

Ce corrects for the difference in event selection efficiency 
(Sec. IIV|) between generic BB events and charmo- 
nium events. It is equal to the efficiency for generic 
BB events divided by the efficiency for the relevant 
charmonium final state. 

Table ^ summarizes the meson yields and efficiencies. 
It also presents the branching fraction product B ■ Be, an 
experimental quantity that does not depend on the sec- 
ondary charmonium branching fractions. There is a 3.1% 
systematic error common to both modes — and, in fact, to 
all final states we study — due to acceptance (1.2%), track 
quality selection (2.4%), uncertainty on Ce (1.1%), and 
number of r(4S') (0.8%). 

The separate e^e^ and fi^ branching fraction mea- 
surements are averaged to obtain the final result (Ta- 
ble^. Each measurement is weighted in the average by 
the inverse of the square of the statistical error plus the 
square of the systematic errors unique to that mode. The 
common systematic error is the largest component of the 
B{B J/ipX) uncertainty. 



VI. Xc PRODUCTION 

We reconstruct Xci and Xc2 mesons in the J/4' 7 decay 
mode. The J/tp mesons must satisfy the criteria listed in 
Sec. IV Al with the additional requirement that the J/tp 
candidate mass m satisfy 3.05 < m < 3.12GeV/c^ for 
e+e~ decays and 3.07 < m < 3.12 GeV/c^ for ■ An 

estimated (74.0±0.4)% of J/V' ^ e+e' and (91.4±0.3)% 
of J/ip — > fJ,+^~ mesons that satisfy all other criteria sat- 
isfy this additional mass selection. 

Photon candidates are EMC clusters in the angu- 
lar range 0.41 < 6 < 2.409 rad with energy between 
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TABLE L Meson yield in on-resonance (20.3 fb^^) and off-resonance (2.6 fb^^) data, and net yield after continuum subtraction. 
ec and Ce are the meson reconstruction efficiency and the relative event selection efficiency; Sn and St are the systematic errors 
on the meson yield and reconstruction efficiency; Stat is the statistical error. The systematic errors on the branching fraction 
products include components unique to that final state. The total uncertainty values include those listed in the "Common" 
rows. Be are the secondary branching fractions and Sb the uncertainty. Tot is the total systematic error (percentage) on the 
B branching fraction to that final state. 





Meson Yield iV^, 


Efficiencies 


Uncertainties (%) 


B Product 


xlO'' 








Mode 


On 


Off 


Net 




Ce 


Sn 


St 


Stat 


Value 


Stat 


Sys 


Be 


Sb % 


Tot % 
































e e 


16095 ± 242 


23± 15 


15914 ± 268 


0.589 


1.02 


4.1 


2.0 


1.7 


650 


±11 


±30 


0.0593 


1.7 


4.9 




13683 ± 154 


67± 18 


13159 ±210 


0.500 


0.99 


0.7 


1.4 


1.6 


615 


±10 


±10 


0.0588 


1.7 


2.3 


Common 


- 


- 


- 


- 


- 


- 


3.1 


- 


- 


- 


3.1% 


— 


0. 


3.1 
































7/ 1 + — 

J/tp e^e 


512 ± 62 


-2 ±3 


528 ± 67 


0.191 


1.04 


6.4 


3.1 


12.7 


68 


±9 


±5 


0.0593 


1.7 


7.3 


J/ip fi+^~ 


614 ± 72 


3±6 


592 ± 86 


0.201 


1.00 


5.4 


2.0 


14.5 


69 


±10 


±4 


0.0588 


1.7 


6.0 


Common 














4.5 








4.5% 


0.316 


10.1 


11.1 


■\y n i T/in '^Z 

Xc2 — » 






























J/ip — > e+e~ 


168 ± 48 


-5 ±3 


210 ± 54 


0.197 


1.04 


3.9 


9.6 


25.7 


26 


±7 


±3 


0.0593 


1.7 


10.3 


J/ip 


208 ± 54 


4± 5 


174 ± 66 


0.207 


1.00 


10.4 


12.1 


38.0 


20 


±8 


±3 


0.0588 


1.7 


16.0 


Common 














5.3 








5.3% 


0.187 


10.7 


11.9 


iP{2S) £+£- 
































573 ± 52 


-6 ±8 


623 ± 81 


0.594 


1.05 


4.3 


1.9 


12.9 


25.8 


±3.3 


±1.2 


0.0078 








437 ± 44 


5± 10 


400 ± 92 


0.535 


1.01 


3.1 


1.6 


23.0 


17.8 


±4.1 


±0.5 


0.0067 






Common 














3.1 








3.1% 






3.1 


ip{2S) J/iPn+ 






























J/ip — > e^e~ 


474 ± 43 


0±2 


476 ± 45 


0.205 


0.99 


2.7 


2.1 


9.5 


53.9 


±5.1 


±1.8 


0.0593 


1.7 


3.8 


J/ijj fi+fi- 


493 ± 42 


0±3 


496 ± 47 


0.215 


0.98 


2.3 


1.4 


9.5 


53.3 


±5.1 


±1.4 


0.0588 


1.7 


3.2 


Common 














3.5 








3.5% 


0.305 


5.2 


6.3 



0.12 GeV and 1.0 GeV. Hadronic showers are suppressed 
by requiring LAT less than 0.8 and A42 < 0.15, while 
clusters from nearby hadronic showers are suppressed by 
requiring that candidates be at least 9° from all charged 
tracks. 

Most photons satisfying these requirements are pro- 
duced in tt'^ decay. We reject a candidate that, when 
combined with any other photon, produces a mass be- 
tween 0.117 GeV/c2 and 0.147 GeV/c^. The second pho- 
ton must have energy greater than 30MeV and LAT< 
0.8, with no requirement on A^i or distance from charged 
tracks. 

A systematic error due to the photon selection crite- 
ria is obtained by comparing the branching ratio t"*" — > 
h^-K^-K^ jr^ — > /i+TT^ in data to that in simulation, where 
is any charged track. We also vary the minimum 
energy requirement from 0.10 to 0.14 GeV and test alter- 
native tt'^ veto regions. We obtain a systematic error of 
3.1% for the Xc\ and 4.4% for the Xci common to both 
the and [i^ [i" final states. An additional compo- 

nent to the systematic error arises from changes in the 
shape of the background, shown in Fig. |21 affecting the 
fit results. It is specific to each final state and mode, and 
amounts to 2.2% (e+e") and 1.3% for the Xci, 

and 9.2% (e+e") and 11.9% {[i^ yr) for the Xc2- 

The photon is constrained to originate at the J/ip ver- 
tex in the calculation of the Xc four-momentum. We 
require -p* < 1.7GeV/c, a requirement that is satisfied by 
Xci or Xc2 mesons from B decays. 



We determine the number of mesons from a fit to the 
plot of the mass difference between the candidate and 
the daughter J/ip masses (Fig. O. We use different sig- 
nal pdfs for the Xci and Xc2- These are formed by con- 
volving the pdf calculated by simulation with a Gaussian 
distribution, where the offset and sigma are constrained 
to be the same for the Xci and Xc2- The background is 
described by a third-order Chebychev polynomial. Sys- 
tematic errors on the fit are obtained as for the J/ip ■ The 
correlation coefficient between the number of Xci and Xc2 
mesons obtained from the fit is 0.19. 

Equation^ is used to determine inclusive B —^ Xci^ 
and B Xc2X branching fractions separately for 
J/ip e+e~ and J/ip — » M^M • In this case. Be = 
l3{Xcj J/-ipl) ■ B{J/ip i^i^), where XcJ is Xci or 
Xc2 and ^ is e or Table summarizes the yields, ef- 
ficiencies, uncertainties and branching fraction products 
BiB XcjX) ■ Bixcj ^ J/V'7) • B{J/iP ^ £+£-). The 
4.5% (xci) and 5.3% (xc2) systematic errors common to 
both e+e^ and lJi~ include photon reconstruction in 
addition to the J/ip reconstruction items. 

As with the J/ip , the inclusive B branching fractions 
to the Xci and Xc2 are calculated separately using the 
J/ip e+e~ and J/ip — > ^J.'^n~ decays. The two values 
are then combined, distinguishing uncertainties common 
to both from those unique to a single final state. 

The branching fraction obtained for the Xc2 is compa- 
rable to that for the Xci and is summarized in Table ITU 
This result is consistent with a prediction from a color 
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FIG. 3: Xci and Xc2 candidates reconstructed in the J/;/' 7 final 
state. Mass difference between the J/i/'7 and J/i/; candidates 
when the J/ip is reconstructed in the (a) e^e~ and (b) 
final states. 



FIG. 4: Mass distribution of ip{2S) candidates reconstructed 
in the (a) e~ and (b) At^^i" final states. 



octet calculation [IQ, and is in contrast to the expecta- 
tion of a null result in a factorization calculation llTj. 



VII. iP{2S) PRODUCTION 

The reconstruction of the tp{2S) in the final 
state is very similar to the J/tp reconstruction outlined 
in Sec. IV Al with the p* requirement tightened to p* < 
1.6GeV/c. Figure 0] shows the resulting candidate mass 
distribution. A fit to extract the number of mesons in 
each plot is performed as for the J/ip , but with the res- 
olution and bremsstrahlung parameters fixed to the val- 
ues found in the higher-statistics J/tp channels. These 
parameters are varied according to their uncertainties as 
one contribution to the systematic error on the fit; the 
remaining contributions are determined as for the J/^ . 

These data are used to calculate the branching frac- 
tion product B{B ^ i^{2S)X) ■ B{tP{2S) and 
are later used in the determination of the p* distribution 
of i/'(2<S') mesons produced in B decay. However, the ex- 
traction of the B 'ijj(2S)X branching fraction requires 
the use of tp{2S) e+e" and ip{2S) — > ^"'"/i" branching 
fractions. Since this same data set has previously been 



used to measure these branching fractions 18l|, we do not 
use ipi2S) -> e+e- events to find B{B -> ^j{2S)X). 

Instead, we use ip{2S) —> J/^tt^tt^ for this purpose. 
The reconstruction oiaip{2S) candidate in this final state 
starts with a J/^^: candidate satisfying the tighter mass 
constraints used in Xcj reconstruction. All charged parti- 
cles, including those failing the "good-track" criteria, are 
assumed to be pion candidates. The pion pair is required 
to be oppositely charged and to have a mass, calculated 
by four- vector addition, in the range 0.45 to 0.60 GeV/c^. 
The mass distribution from simulation is compared to the 
measured 19] distribution to obtain a systematic error of 
0.5% on reconstruction efficiency. Finally, the p* of the 
'0(25') candidate is required to be less than 1.6GeV/c. 

Figure |S1 displays the mass difference between the 
^{2S) and the J/ip candidates separately for J/^ e+e~ 
and J/ip — > /J."*"^". As for the other final states, the 
distributions are fit to obtain the number of mesons. 
The resolution smearing parameters are not required 
to be the same for the two plots, but are consistent: 
1.5 ± 0.8MeV/c2 (e+e") and 1.8 ± 0.5 MeV/c^ (/x+Ai"). 
The secondary branching fractions in Eq. ^ are in this 
case Be = 6(^(25') ^ J/^pn+Tr-) ■ B{J/ip £+£"). 
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FIG. 5: tp{2S) candidates reconstructed in the J/iI^-k^-k' final 
state. Mass difference between tlie J/ip-R^-K~ and J/ip candi- 
dates when the J/^ is reconstructed in the (a) e^e~ and (b) 
final states. 



VIII. DIRECT BRANCHING FRACTIONS 

To obtain the branching fraction for J/ip mesons pro- 
duced directly in the decay of B mesons, we subtract the 
feeddown contributions to the inclusive branching frac- 
tion due to the decay of Xci, Xc2, and ^p(2S) mesons. 
For the Xci and Xc2i the feeddown branching fraction is 
B{B -> XcjX) ■ Bixcj ^ J/i'l), while for the 7/;(25), it 
is B{B ^ ^p{2S)X) ■ B{^p{2S) J/i^X). 

Similarly, the feeddown from the il>{2S) to the Xci and 
Xc2 is B{B ^ ^{2S)X) ■ B{^{2S) ^ XcJl)- 

Note that a number of uncertainties are common to 
both the inclusive and feeddown components, including 
track quality and particle identification criteria, B{xci — > 
J/V't), and B{xc2 J/i^l)- We use world average val- 
ues [13 for the '0(25') branching fractions. The resulting 
direct branching fractions are summarized in Table Hll 



IX. p* DISTRIBUTIONS 

The momentum distributions of charmonium mesons 
provide an insight into their production mechanisms. 
Since we do not fully reconstruct the B meson, we can- 



TABLE II: Summary of B branching fractions (percent) to 
charmonium mesons with statistical and systematic uncer- 
tainties. The direct branching fraction is also listed, where 
appropriate. Last column contains the world average values 



Meson 


Value 


Stat 


Sys 


World Average 




1.057 


±0.012 


±0.040 


1.15 ±0.06 


J/ip direct 0.740 ±0.023 ±0.043 


0.80 ± 0.08 


Xci 


0.367 


±0.035 


±0.044 


0.36 ± 0.05 


Xci direct 


0.341 


±0.035 


±0.042 


0.33 ± 0.05 


Xc2 


0.210 


±0.045 


±0.031 


0.07 ±0.04 


Xc2 direct 


0.190 


±0.045 


±0.029 




V(25) 


0.297 


±0.020 


±0.020 


0.35 ±0.05 



not determine the meson momentum in the B rest frame 
and instead use p*, the value in the T{AS) center-of-mass 
frame. The difference, due to the motion of the B in the 
center-of-mass frame, has an rms spread of 0.12 GeV/c. 



A. Inclusive p* Distributions 

To measure the p* distributions of J/V', Xci, Xc2, and 
ip{2S) mesons produced in B decays, we create mass 
or mass-difference histograms of on-resonance candidates 
with p* in the desired range. The e~^e~ and jX^ n~ final 
states are again treated separately. The distributions are 
then fit, with all signal pdf parameters (other than the 
number of mesons) fixed to the values obtained from the 
earlier fits. The fits are performed for lOOMeV/c wide 
p* ranges, and in each case the sum of the yields differs 
from the original fit by fewer than ten events. 

In the case of the J/tp , we perform similar fits on the 
off-resonance data and perform a continuum subtraction 
for each p* bin. Since there are no statistically signifi- 
cant off-resonance Xci, Xc2, or ip{2S) signals, we do not 
perform a continuum subtraction in these cases. 

The yield in each bin is corrected by the reconstruction 
efficiency obtained from simulated data, which decreases 
by approximately 10% between and 2 GeV/c. The yield 
is then multiplied by an overall normalization factor for 
that particular final state and mode, which adjusts the 
sum of all bins to the earlier branching fraction measure- 
ment. We then perform a weighted average of the two 
distributions for the J/ip, Xci) or Xc2, or the four distri- 
butions for the ip{2S), to obtain the distributions shown 
in Fig.EHSl For this purpose, we use the ^{2S) e+e" 
and i^{2S) fJ^^fJ-^ branching fractions from Ref. [l^ . 
In all cases, the distributions that are combined are con- 
sistent within statistical errors. 
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FIG. 6: B J/^X branching fraction as a function of p* 
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FIG. 8: B — > tp{2S)X branching fraction as a function of p* 
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FIG. 7: Branching fraction as a function of p* for (a) B 
XciX and (b) B X^^X. The distribution includes a small 
feeddown component from the ip{2S) (solid curve). 



B. Direct p* Distributions 

The J/ip p* distribution (Fig. |SJ) includes components 
due to mesons from the decays Xci J/i'l, Xc2 J/i'l, 
and ^{2S) — * J/tpX. To measure these distributions, we 
repeat the analysis with the data binned by the p* of the 
J/V' daughter. The resulting J/^p feeddown distributions 
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FIG. 9: Contributions to the B — » J/tpX branching fraction 
as a function of p* due to feeddown from (a) Xci, (b) Xc2 
and (c) ip{2S) mesons. 



are presented in Fig. |3 

Note that we are using only the J/?/;7r+7r^ decay mode 
to obtain the J/'ip distribution from ip{2S) decay. In fact, 
10.5% of ip{2S) J/ipX decays are modes other than 
J/ipnn. If we instead use the simulated J/ip distribution 
for this 10.5%, Fig. changes by no more than a small 
fraction of the statistical error bar in any bin. 

Subtracting these three components from the inclusive 
J/ip distribution in Fig. |H| leaves the contribution due to 
the J/4> mesons produced directly in B decay (Fig. I10|). 

The superimposed histogram is a calculation of the ex- 
pected distribution, which includes color octet and color 
singlet components. We use a recent NRQCD calcula- 
tion for the color octet component. The authors 
attribute the singlet component to J/ipK^*^ production, 
which we obtain from simulation. The two are normal- 
ized to obtain the best fit to our data. Possible sources 
of the apparent excess at low momentum are an intrinsic 
charm component of the B |21l , the production, together 
with the J/'ip , of baryons |22l| , or an sdg hybrid [23j . 

The small feeddown contribution to Xci and Xc2 from 
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(points) . The histogram is the sum of the color-octet compo- 
nent from a recent NRQCD calculation 20] (dashed line) and 
the color-singlet J/t/jJ^^*^ component from simulation (dotted 
line). 
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FIG. 12: Helicity distribution of J/ip mesons produced in the 
decay of (a) Xci, (b) Xc2, and (c) tp{2S) mesons. 
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FIG. 11: Helicity of J/tp mesons produced in B decay with 
p* > 1.1 GeV/c (dots) and p* < 1.1 GeV/c (open squares). 



'ip{2S) decay is calculated by simulation and is shown in 
Fig.0 



We proceed as for the J/ip p* distribution, with data 
categorized into ranges of width 0.1 in cos6h for two 
different momentum ranges, which we choose as p* < 
1.1 GeV/c and 1.1 < p* < 2.0GeV/c. We fit the on- and 
off-resonance mass distributions to obtain yields in each 
bin and perform a continuum subtraction. We correct 
using the reconstruction efficiency obtained from simula- 
tion for that range, although we observe little dependence 
of efficiency on helicity. We then apply separate normal- 
ization factors to the e~^e~ and fi^ data such that the 
total branching fraction (summed over the two p* ranges) 
agrees with the value obtained earlier for that mode. The 
distributions from e'^e~ and are consistent and are 

averaged to obtain the helicity distributions for each of 
the two p* ranges (Fig. Illfl. 

We fit each distribution with a function 1 -I- a • cos^ Oh 
to obtain the polarization a, where a = indicates the 
sample is unpolarized, a = 1 transversely polarized, and 
a = — 1 longitudinally polarized. The high p* region, 
which includes the two-body B decays, is more highly 
polarized, a — —0.592 ± 0.032, than the lower p* region, 
a = -0.196 ±0.044. 

We assign a systematic error of 0.008 to these polariza- 
tions by instead considering the reconstruction efficiency 
to be independent of helicity. 



X. J/V- HELICITY 

The helicity 6h of a J/i/i candidate is the 

angle, measured in the J/tp rest frame, between the pos- 
itively charged lepton and the flight direction of the J/ip 
in the T(4S') center-of-mass frame. A more natural def- 
inition would use the B rest frame, but it cannot be de- 
termined in this analysis. Simulation indicates that the 
rms spread of the difference between the two definitions 
is 0.085 in cos 9h- 



B. Direct J/tp Helicity 

We determine the helicity distributions of J/^ mesons 
produced in the decay of Xci, Xc2, and ?/'(2S') in the same 
way we calculate the p* feeddown. Because of the limited 
statistics of these samples, we combine the two momen- 
tum regions used in the inclusive analysis. The resulting 
feeddown helicity distributions are shown together with 
the polarization fits in Fig. ^1 We subtract these from 
the sum of the two distributions in Fig. ^] to obtain the 
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FIG. 13: Helicity distribution of J/i/; mesons produced di- 
rectly in the decay of B mesons. 



helicity distribution for the J/ip produced directly in B 
decay (Fig. [TSJl. The polarization, a = -0.46 ± 0.06, is 
slightly out of the range —0.33 to 0.05 predicted by an 
NRQCD calculation [2j, but other authors have argued 
[2 fx] that relativistic corrections reduce the reliability of 
the calculation. The systematic uncertainty of 0.008 ob- 
tained above is small compared to the statistical error. 
This result is difficult to compare directly with that from 
CDF 2.6|, due to the different mixture of b mesons and 
baryons, and the distinction between the effective helicity 
calculated there and the true helicity. 



XI. SUMMARY 

We have reported new measurements of B meson de- 
cays to final states including charmonium mesons, which 
are summarized in Table ITTl We have presented a num- 
ber of momentum distributions. The distributions of the 
feeddown J/^p daughters of the Xci, Xc2, and iIj{2S) have 
not previously been measured, and allow us to more ac- 
curately determine the distribution for J/^ mesons pro- 
duced directly in B decay. The direct J/ip distribution 
is compared to a recent NRQCD calculation and appears 
to indicate an excess at low momentum. 

The J/iJj helicity distribution, which has also has not 
previously been published, indicates that the polarization 
of direct J/ip mesons is slightly out of the range predicted 
by an NRQCD calculation. 
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